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Spatial Point Analysis of Quantum Dot Nucleation Sites on InAs Wetting Layer 
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We perform spatial point analysis of InAs quantum dot nucleation sites and surface reconstruction 
domain pattern on an InAs wetting layer, giving insights for quantum dot nucleation mechanism. An 
InAs wetting layer grown to 1.5 monolayers in thickness on a GaAs(OOl) substrate has been observed 
at 300 °C by using in situ scanning tunneling microscopy. The surface exhibits (1 x 3)/(2 x 3) and 
(2x4) reconstruction domains. A nearest-neighbor analysis finds that point pattern of quantum dot 
precursors was more similar to that of (1 x 3)/(2 x 3) domains which are specific to Ga-rich region. 
This provides the evidence that InAs quantum dot nucleation is induced by Ga-rich fluctuation 
within an InAs wetting layer. 

PACS numbers: 68.37.Ef, 68.43.Hn, 68.47.Fg, 68.55.ag 
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Quantum dots (QDs) are potentially used in semicon- 
ductor laser devices and single photon sources of quan- 
tum computation and quantum communication [l|-[4J]. 
Although it has been pointed out that highly dense and 
uniform QD arrays are essential for the efficiency of the 
devices, little is known of the growth mechanism of QDs 
to control the nucleation sites on a Stranski-Krastanow 
(SK) grown wetting layer (WL). Some atomic-level the- 
oretical studies on dynamics of surface atoms have been 
carried out to understand the growth mechanisms [5|- 
Q. First principle calculations showed that the migra- 
tion barrier energy of In adatom on GaAs(OOl) surface is 
higher than that on lML-InAs/GaAs(001) [1,0. Using 
kinetic Monte Carlo (kMC) simulations [8], Tsukamoto et 
at. found that some migrating In adatoms were captured 
on Ga-rich fluctuation, within an In/Ga mixed layer, to 
become a nucleation site [9] . To the best of our knowl- 
edge, however, there has not been reported any direct 
evidence that alloy fluctuation becomes a QD nucleation 
site. It is still vital to investigate WL surface in an atomic 
scale, in particular the surface reconstruction, preceding 
QD nucleation. 

Since surface reconstruction changes microscopically 
and dynamically in the course of WL growth J10l-[l2| , in 
situ scanning tunneling microscopy (STM) during molec- 
ular beam epitaxy (MBE) growth at high temperatures, 
such as STMBE [13], is one of powerful tools to ob- 
serve it. It is reported that fast Fourier transform anal- 
ysis of atomic-scale in situ STM images of InAs WL 
on a GaAs(OOl) substrate, as well as reflectance high- 
energy electron diffraction (RHEED) measurements, has 
revealed that the surface reconstruction changes from 
c(4 x 4) to the mixed structure of (1 x 3)/(2 x 3) do- 
mains and (2 x 4) domains prior to QD formation [9]. 

Structure models of (1 x 3)/(2 x 3) and (2 x 
4) surface reconstructions have been investigated by 
many researchers using core-level photoemission spec- 
troscopy [14]], reflectance-difference spectroscopy [l5j, Q>b 
initio calculations in a local density approximation p, 
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FIG. 1. (Color online) Schematic diagrams of surface recon- 
struction models of InAs/GaAs wetting layer on GaAs(OOl) 

reproduced from literature: (a) (1 x 3) [lj], (b) (1 x 3) flak 

(c) (2 x 3) 6], (d) a2(2 x 4) EElH, (e) £2(2 x 4) [IE El 
Background shade indicates unit cell. 



LL5[, and STM observations [l6|, [l3l 5 which are, however, 
still under discussion. Figure [1] shows the schematic dia- 
grams, reproduced from the literature, of some represen- 
tative surface reconstruction models, in which only the 
atoms near the surface are illustrated. The unit cells of 
(1x3) and (2 x 3) have one or two of Ga atoms near 
the surface whereas those of (2 x 4) have none. In other 
words, Ga-rich fluctuation in InAs/GaAs WL, which is 
expected to become QD nucleation sites, likely forms 
(1 x 3)/ (2 x 3) surface reconstruction domains. It is 
crucial to investigate the relationship between QD nu- 
cleation sites and surface reconstruction domains, but no 
researcher has done yet. It is quite challenging, however, 
to witness QD formation specific to a certain reconstruc- 
tion domains using typical STM scanning speed because 
QD formation occurs on dynamically changing recon- 
struction, over a timescale of a few seconds even at very 
low WL growth rates. In this paper, we firstly demon- 
strate a statistical approach to this problem, namely the 
statistical comparison of the distribution of surface re- 
construction domains and that of QD nucleation sites. 

The distribution of reconstruction domains and QD 
nucleation sites is characterized by spatial point patterns; 



that is a regular (ordered) pattern, a Poisson (random) 
pattern, and a clustered (aggregated) pattern [18|, fl9J . 
In a regular pattern, points are distributed uniformly. A 
Voronoi tessellation, that is polygonal decomposition of a 
space by perpendicular bisector lines among neighboring 
points, is often used in spatial point analysis. The stan- 
dard deviation of Voronoi cell area, <jv c , represents point 
patterns. Here, we have to be careful of the Voronoi cells 
touching the boundary of the study region because of 
"edge effect" that contribution from points outside can- 
not be taken into account. In this study, such invalid 
Voronoi cells were excluded from the study region. For 
more precise analysis, second-order properties of point 
patterns like nearest-neighbor distances are useful J19J - 
[22[. Let r denote the distance to the nearest point from 
a randomly selected location in the study region R. The 
F(d) function denotes the cumulative frequency distribu- 
tion of r [19j and hence the probability that r occurs less 
than any particular distance d. Since the F function is 
practically identical to the probability of plotting a ran- 
dom point within any of circles Ci(d), of each radius d, 
centered on each of the points, it is simply computed by 



F(d) 



Area(U i Ci(rf)n J R) 
Area (R) '' 



(1) 



where the numerator is the area of the union of the cir- 
cles, and the denominator is the area of the study region 
R. To compare F(d) between different study regions, d 
should be normalized by the factor / as follows: 



f = y/Aie&(R)/N, 



(2) 



where N is the number of points [20] . 

A piece (11 x 13 x 0.6 mm 3 ) of GaAs(OOl) crystal was 
used as a substrate. First the surface was thermally 
cleaned to remove the oxide layer under 1 x 10 -4 Pa of 
an AS4 atmosphere in an MBE growth chamber. Next, 
a GaAs buffer layer was grown on the surface by using 
MBE until atomically smooth surface was obtained. The 
substrate was annealed at 430 °C for 0.5 h to confirm the 
formation of c(4 x 4) reconstruction with RHEED. An 
STM unit was transferred to the sample holder in the 
growth chamber and started scanning. This means that 
the sample was neither cooled to room temperature nor 
exposed to air to be observed with STM. The tip bias was 
-3 V and the tunnel current was 0.2 nA. A flux of In was 
irradiated at the InAs growth rate of 2.5 x 10 -4 MLs -1 . 
After 1.5 monolayer (ML) of InAs WL growth, the sub- 
strate temperature was decreased to 300 °C and the AS4 
flux was shut off so that the surface reconstruction should 
be observed in an atomic scale. Another samples was pre- 
pared in the same way but SK growth of InAs WL was 
continued at 430 °C until QD precursors were formed [9]. 
An STM image of QD precursors was recorded to analyze 
the distribution of QD nucleation sites. 

Figure [2] shows the filled-state STM image of 1.5 ML 
of InAs WL recorded at 300 °C. The image seems more 
or less scratched because of migrating In adatoms on the 
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FIG. 2. (Color) 30 nm x 30 nm filled-state STM image of 
InAs WL on GaAs(OOl). 
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FIG. 3. (Color) Map of As-dimer stripe pitch along [110] 
azimuth measured in each cell of 25 x 25 mesh in Fig. [2] 
Schematic diagram of magnified mesh cell is also illustrated 
aside to show the relationship between stripe pitch measured 
and corresponding unit cell of InAs WL surface reconstruction 
in Fig. [U 



surface but shows the stripes due to As dimers clearly 
enough to measure the pitch. The stripe pitch corre- 
sponds to the unit cell length along [110] azimuth of InAs 
WL surface reconstructions in Fig. [TJ Although it is dif- 
ficult to discuss the unit cell size of strained and mixed 
surface reconstructions, it is expected to have some in- 
termediate value between those of GaAs and InAs. We 
assumed that the stripe pitch along [110] azimuth ranged 
0.6-1.0 nm for c(4 x 4), 1.0-1.4 nm for (1 x 3)/(2 x 3), 
and 1.4-2.0 nm for (2x4). Before measuring the stripe 
pitch, the STM image [Fig. [2] was divided by a 25 x 25 
mesh. As shown in the schematic diagram in Fig. [3j the 
size of each mesh cell is 1.2 nm which is comparable to 
the unit cell sizes of the InAs WL surface reconstruc- 
tions. The stripe pitch was measured from the STM line 





(a) (1x3)/ (2x3) domains 
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(b) (2x4) domains 

FIG. 4. (Color) Surface reconstruction domains, indicated by 
oval markers, and Voronoi tessellations of (a) (1 x 3)/(2 x 3) 
and (b) (2x4) surface reconstruction maps. 



profile along [110] azimuth for each mesh cell. The data 
are plotted in the color map of Fig. 03 Most cells show 
(1 x 3)/ (2 x 3) or (2x4) surface reconstruction although 
some cells are blank because of step edges. 

It is predicted that Ga-rich fluctuation, comprised of at 
least eight Ga atoms in an In/Ga mixing layer, should be 
a QD nucleation site according to kMC simulations |3|. 
For such Ga-rich fluctuation to be formed, four unit cells 
of (1 x 3)/ (2 x 3) need to be contiguous, each of which 
has one or two Ga atoms. In this study, groups of four 
neighboring mesh cells, having the same surface recon- 
struction, were located in the map [Fig. [3] and indicated 
by oval markers in Fig. |U For each of(lx3)/(2x3) and 
(2 x 4) surface reconstructions, the centroid points of the 
domains were marked and their coordinates were mea- 
sured by using Image J software J23, [24|- The centroid 
coordinates were used for the Voronoi tessellations of the 
reconstruction domain maps [Fig. 0] and computation of 
the F function [13. 

Figure [5] shows the 150 nm x 150 nm STM image 
of InAs QD precursors immediately after nucleation [9| . 
The STM image was divided into 3x3 regions. For 
each region, measurement of QD coordinates, a Voronoi 
tessellation, and computation of the F function were per- 
formed as well. 

Table. U lists the density and standard deviation of 
valid Voronoi cells computed for the surface reconstruc- 
tion domains in Fig. |4] and the QD precursors in Fig. [5] 
Each standard deviation is normalized by each study re- 
gion area. For QD precursors, the minimum and maxi- 




FIG. 5. (Color) 150 nm x 150 nm STM image of QD pre- 
cursors nucleated on SK grown InAs WL at 430 °C [9| and 
Voronoi tessellations computed for each of 3 x 3 regions. 



TABLE I. Density p and standard deviation of Voronoi cell 
area ay c of surface reconstruction domains and QD precur- 
sors. 



pqo 1 



2 ) 



OVc 



(1 x 3)/(2 x 3) domains 6.2 0.38 

(2 x 4) domains 4.2 0.31 

QD precursors 0.96-1.7 0.20-0.59 



mum data in the nine regions are shown. The densities 
of surface reconstruction domains were similar to those 
of QD precursors. The standard deviations of surface 
reconstruction domains were in the range of QD precur- 
sors. The similarity in these properties implies some re- 
lationship between surface reconstruction domains and 
QD precursors. 

The F function will give more precise information. Fig- 
ure [6] shows the F(d) traces calculated for the surface 
reconstruction domains [H] and the F(d) envelope of the 
QD precursors [[5]. The F(d) envelope of typical Poisson 
patterns was calculated by accumulating 50 simulations 
of scattering 50 random points. 

Both traces of (1 x 3)/ (2 x 3) and (2 x 4) domains 
were similar to the envelope of QD precursors although 
they differs in the detail. The traces of the surface re- 
construction domains and the QD precursors were plot- 
ted between those of the ordered pattern and the Pois- 
son patterns. This shows that surface reconstruction do- 
mains and QD precursors are distributed in a rather or- 
dered pattern than a random pattern. An ordered pat- 
tern likely occurs when there is repulsive force among 
points. It is difficult to discuss the origin of repulsive 
force by STM images, one can consider the surface strain 
distribution. First principles calculations of WL would 
give some insights for the surface stain distribution. It 
is possible that the uniformity of QD nucleation, which 
is rather ordered, is originated from the distribution of 
surface reconstruction domains. 

The trace of (2 x 4) domains is located rather ordered 
in the QD envelope, and completely deviates from the 
QD envelope in small d region as can be seen in the mag- 
nified view in Fig. [6j The trend of F(d) where d is small 
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FIG. 6. (Color) Traces of F functions of surface reconstruc- 
tion domains and typical ordered point pattern. Envelope 
regions of QDs in nine regions and typical Poisson patterns 
by accumulating 50 simulations are also shown. Magnified 
view of region where d is 0.36 to 0.40 is superimposed. 



is dependent on the pattern in dense-point areas. On the 
other hand, the trend of F(d) where d is large is depen- 
dent on the pattern in sparse-point areas. The pattern in 
dense-point areas is particularly crucial because of high 
surface strain which severely affects the distribution. In 
the sense of not deviating in dense-point areas, the trace 



of (1 x 3)/ (2 x 3) domains represents the QD precur- 
sor envelope better than that of (2 x 4) domains. In 
other words, the distribution of QD precursors are more 
similar to that of (1 x 3)/ (2 x 3) domains than that of 
(2 x 4) domains. The similarity in the density and the 
point pattern implies that QD nucleation is related to 
(1 x 3)/ (2 x 3) domains where Ga-rich fluctuation is as- 
sumed to be formed. This also means that QD nucleation 
sites are already determined at the time when InAs WL 
is grown to 1.5 ML or before. 

In conclusion, we have shown, by using in situ STM, 
the similarity in the density and the spatial point pat- 
terns between (1 x 3)/ (2 x 3) surface reconstruction do- 
mains on 1.5 ML of SK grown InAs WL and QD precur- 
sors nucleated. This implies that QD nucleation site is 
related to the distribution of (1 x 3)/ (2 x 3) domains at the 
stage of 1.5 ML-grown InAs WL. Since a (1 x 3)/(2 x 3) 
domain composed of four unit cells is assumed to con- 
tain Ga-rich fluctuation with at least eight Ga atoms, 
this model provides a consistent evidence that QD nucle- 
ation is induced by tiny alloy fluctuation as predicted by 
kMC simulations [9] . The mechanism of QD nucleation, 
exhibited here, has important technological implications 
for the self-assembly and the artificial arrangement of 
QDs. 
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